Tunnel magnetoresistance of quantum dots coupled to ferromagnetic leads 
in the sequential and cotunneling regimes 
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We study electronic transport through quantum dots weakly coupled to ferromagnetic leads with 
collinear magnetization directions. Tunneling contributions of first and second order in the tunnel- 
coupling strength are taken into account. We analyze the tunnel magnetoresistance (TMR) for all 
combinations of linear and nonlinear response, at or off resonance, with an even or odd dot-electron 
number. Different mechanisms for transport and spin accumulation the various regimes give rise to 
different TMR behavior. 
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I. INTRODUCTION 

The study of spin-polarized electron transport through 
nanostructures with strong Coulomb interaction is a rela- 
tively new field of theoretical and experimental research, 
residingin the intersection of the fields of spintronics 
[B0E1 H and transport through nanostructures 0,IS0> 
respectively. The interplay of finite spin polarization and 
Coulomb blockade gives rise to a complex transport be- 
havior in which both the electrons' charge and spin de- 
gree of freedom play a role A convenient minimal 
model system to study this interplay consists of a single- 
level quantum dot coupled through tunnel barriers to 
ferromagnetic electrodes. Experimentally such systems 
may be realized in various ways, including self-assembled 
dots in ferromagnetic semiconductors [!j , ultrasmall alu- 
minum nanoparticlesjnj , carbon nanotubes 0, 0, ^| , 
or single molecules |l4j . 

The properties of spin-polarized transport through sin- 
gle magnetic tunnel junctions have already proven tech- 
nological relevance in information-storage devices based 
on the tunnel magnetoresistance (TMR) effect, i.e., the 
observation that the current flowing through the junction 
depends on the relative orientation of the leads' magne- 
tizations. It is maximal for the parallel and minimal for 
the antiparallel configuration. Quantitatively, it can be 
characterized by 



TMR 



Ip — Iai 
Iap 



(1) 



where Ip and Iap are the currents for the parallel and an- 
tiparallel configuration, respectively. Julliere found 0] 
that the TMR for a single tunnel junction is related 
to the degree p of spin polarization of the leads' den- 
sity of states, p={p+ - p~)/(p + + p~), by TMR Ju11 = 
2p 2 /(l — p 2 ), where p + and p~ are the spin-majority and 
spin-minority densities of states in the electrodes, respec- 
tively. Julliere's formula immediately follows from the 
fact that the transmission probability of an electron with 



spin a through the barrier is proportional to the product 
of the (spin-dependent) densities of states for spin a in 
source and drain. 

Once a nanoscopic island is placed in between the fer- 
romagnetic leads the situation becomes much more com- 
plex for two reasons. First, there are different types of 
transport processes that depend on the leads' spin po- 
larization in a different manner, such as sequential tun- 
neling, non-spin- flip, and spin-flip cotunneling (for non- 
spin-flip cotunneling an electron of given spin is trans- 
ferred through the system, while for spin-flip cotunnel- 
ing both the spin of the transferred electron as well as 
the dot spin changes during the process). Second, a 
non-equilibrium spin accumulation can partially polarize 
the island, which, in turn, affects the total transmission 
through the device. Therefore, the TMR will, in general, 
deviate from Julliere's value. It will, furthermore, be dif- 
ferent for different transport regimes. The measurement 
of the TMR as a function of temperature, bias and gate 
voltages, will, thus, reveal information about the under- 
lying transport processes as well as the spin accumulation 
on the island. 

Spin-dependent transport through a single-level quan- 
tum dot in the sequential-tunneling regime with 
collinearly magnetized leads has been analyzed in 
Refs.IHIllH This has been extended [IE EO, Ell E3 to 
noncollinear configurations with arbitrary relative angle, 
for which a precession of the dot spin about an intrinsic 
exchange field gives rise to non-trivial dependence of the 
angle-dependent conductance. In the present paper, we 
analyze the TMR for collinear magnetization beyond se- 
quential tunneling. This covers the Coulomb-blockade 
regime, in which sequential tunneling is exponentially 
suppressed, and transport is dominated by cotunneling 
El IM Hi HI |H |2l Hj. But even when sequen- 
tial tunneling is possible, second-order corrections to the 
current become important for increasing tunnel- coupling 
strengths. This includes the above-mentioned cotunnel- 
ing processes but also terms associated with renormaliza- 
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FIG. 1: Single-level quantum dot coupled to ferromagnetic 
leads. The magnetic moments of the electrodes are either 
parallel or antiparallel to each other. 

tion of level position and tunnel-coupling strength |30| . 
Recently, we studied spin-dependent transport for a spe- 
cific transport regime, namely, cotunneling deep inside 
the Coulomb-blockade valley [3l| - 

Our objective for the present paper is to analyze the 
TMR in the full parameter space defined by the gate and 
bias voltages. This includes the linear- and nonlinear- 
response regime as well as the cases of even and odd 
dot occupation. We find that the TMR reaches Jul- 
liere's value only when the transport is fully carried by 
non-spin-flip cotunneling. This happens in the Coulomb- 
blockade valleys in which the dot is either empty or dou- 
bly occupied, where the dot remains unpolarized, as well 
as for large bias voltage in the Coulomb-blockade valley 
with an odd dot-electron number. For all other regimes, 
though, the TMR is reduced below Julliere's value. 

II. MODEL 

We consider transport through a single-level quantum 
dot. The dot is coupled to two ferromagnetic electrodes 
with collinear, i.e., either parallel or antiparallel, magne- 
tizations, see Fig. ^ The dot level e can be tuned by 
a gate voltage, but is independent of the symmetrically- 
applied transport voltage. 

We model the system by an Anderson-like Hamiltonian 
of the form 

H = H L + H R + H U + H T . (2) 

The first and second terms represent the left and 
right reservoirs of noninteracting electrons, H r — 
J2 q a £ rqv c lqa c rq<7, for r = L, R, where c\ qa (c rqa ) is the 
creation (annihilation) operator of an electron with wave 
number q and spin a in the lead r, whereas e rqa denotes 
the corresponding single-particle energy. The dot is rep- 
resented by 

H D = £dU<r + Ud\d^d\di, (3) 
<r=T4 

with d a [da) creating (annihilating) an electron on the 
dot with spin a and energy e, and U is the charging 



energy for double occupancy. There are four possible 
states for the quantum dot: empty dot (x = 0), singly- 
occupied dot with a spin-up (\ — T) or spin-down (\ = j) 
electron, and doubly-occupied dot (x — d). Tunneling 
between dot and leads is described by 

Ht = Y Y {t rqa c\ qa d a + t* rqa d a c rqa ) , (4) 
r=L,R go- 
where t rqa are the tunnel matrix elements. Tunnel- 
ing gives rise to an intrinsic broadening T a of the dot 
levels, given by the Fermi-golden-rule expression T a = 
Er^L.R 1 ?. with = 2nY /q \t rq v\ 2 5(u - e rqa ). As- 
suming the matrix elements t rqa to be independent of 
the wave number and spin orientation, we get Y a r = 
2w\t r \ 2 Pr, with p% denoting the spin-dependent density 
of states in lead r. In the following we assume the latter 
to be independent of energy within the electron band. 
Furthermore, we introduce the degree of spin polariza- 
tion p r = (p+ — Pr)/{pt + Pr) °f ^ ea d r ■, an d express 
the four respective couplings in terms of spin polariza- 
tion as r+ (_) = T r (l±p r ), where T r = (r++T-)/2. In 
general, the leads may have different spin polarizations 
and/or coupling strengths to the dot. In the following, 
however, we assume j»l = Pr = V and Tl = Tr = T/2. In 
the weak coupling regime, typical values of the dot-lead 
coupling strength T are of the order of tens of peV (23 . 

III. METHOD AND TRANSPORT EQUATIONS 

We calculate the transport properties of the system 
by ma king use of a real-time diagrammatic technique 
|30ll32ll33l |. Its main idea is to integrate out the electronic 
degrees of freedom in the leads in order to arrive at an 
effective description of the dot subsystem. The dynam- 
ics of the subsystem is then described by a reduced, four- 
dimensional, density matrix with density matrix elements 
P* 1 (t). The time evolution of the reduced system can be 
represented graphically as a sequence of irreducible dia- 
grams on the Keldysh contour. An example of such time 
evolution is shown in Fig. |2 where the upper and lower 
branches of the Keldysh contour represent the forward 
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FIG. 2: An example for the time evolution of the reduced 
density matrix. The grey regions define irreducible diagrams 
of first and second order in tunneling, respectively. The di- 
rection of each tunneling line indicates whether an electron 
of respective spin leaves or enters the dot, thus, leading to 
a change of the dot state, as indicated on the forward and 
backward Keldysh propagators. 
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and backward propagators. Tunneling is represented by 
vertices, that are connected in pairs by tunnel lines. Each 
grey region in Fig. defines an irreducible diagram that 
corresponds to a transition of the dot state. First- and 
second-order transport in the tunnel- coupling strength T 
is described by diagrams containing one or two tunnel 
lines, respectively. Since we consider only collinear mag- 
netic configurations of the leads and tunneling is spin 
conserving, the natural choice of the spin-quantization 
axis results in vanishing of all non-diagonal density ma- 
trix elements, and only the diagonal ones, = P x , need 
to be considered. They are nothing but the probability 
to find the dot in state x- 

The time evolution of the reduced density matrix is 
governed by a generalized master equation |30| that in 
the stationary limit reduces to 



To determine P^ and Px', we have to expand the 



(i) 



XX r X ' 



(5) 



master equation, Eq. (JSJ, order by order, 



(12) 
(13) 



The evaluation of P^ ] and P x 1] from Eqs. and O 
has to be done with some care. As we will see below, we 
have to distinguish between the two cases in which se- 
quential tunneling is either present or exponentially sup- 
pressed. 



Perturbation expansion in the presence of 
sequential tunneling 



where E x ' x describes the irreducible diagram parts with 
transitions from state \ to x' ■ The electric current is 
given by 



le 
2h 



xx' 



p 

x r x 1 



(6) 



where the self energy SL is modified as compared to 



J x'x 



to account for the number of electrons transferred 



through the barriers. The rules to calculate £ X ' X and 
SL X are given in the appendix. 

Our goal is to calculate the current up to second or- 
der in the tunnel- coupling strength T. For this, we first 
expand the self-energies £ X ' X and £ x / x order by order, 



J x'x 



(7) 



where the order corresponds to the number of tunnel lines 
of a diagram. Consequently, the entire problem is re- 
duced to the calculation of all the self-energies with the 
aid of the diagrammatic rules. 

For an accurate perturbation expansion of the current, 
we also need to expand the probabilities in orders of T, 



p(0) 

X 



p(l) 

X 



with the normalization condition 



(m) 



S m ,0 



(8) 



(9) 



The first- and second-order contributions to the current 
are then given by 



In regime where the sequential tunneling is allowed, 
one can use the perturbation expansion presented in the 
previous subsection. In particular, one can determine the 
zeroth-order probabilities P x from Eq. H12|) and, then, 
plug the result into Eq. (|13|) in order to evaluate the 
first-order corrections P x ^ . Having calculated the prob- 
abilities, one can use the result to get the current from 
Eqs. (|10|l and Ijllfl in first and second order, respectively. 



B. Perturbation expansion in the 
Coulomb-blockade regime 

In the Coulomb-blockade regime, several of the first- 
order self-energies are exponentially small as they 
are associated with energetically forbidden sequential- 
tunneling rates. As a consequence, all addends in the 
first-order master equation, Eq. I|12|) . are exponentially 

small: either the state x is classically forbidden, i.e., P x ^ 
is exponentially suppressed, or the state x is classically 
allowed but then the corresponding self energies Ey/L are 
exponentially small. 

This is not a problem for the Coulomb-blockade valleys 
with an even number of electrons, k^T, \eV\ <C s, £ + U 
and kftT, \eV\ <C —e,—e — U, since for this case, the 
first-order master equation, Eq. ijT^t , yields P^ = <5 x ,o 
and P x ^ = <5 Xi d, respectively, i.e., there is only one 
classically-allowed dot state. The situation is different 
for the Coulomb-blockade valley with an odd number of 
electrons, k^T, \eV\ <§; — e, e + U, where both x = J \ an d 
X = i are classically occupied. In this case, Eq. i(T2l 
simplifies to 



/W = 

j(2) = 



xx' 



2h ^ L x 



2)p(0) , yJ(l)p(l) 
X X ^ x'x X 



(10) 

(11) 



xx' 



(1) 

no 

M 







.-(1) 



/ ^o (0) \ 



(0) 
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(0) 

1 

(0) 



(14) 
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i.e., we obtain P. 



(0) 



,(0) 



while the individual 



occupations p| ' and Pj ^ remain undetermined. Fur- 
thermore, we find that P^ and P^ 1 ' drop out of the 
second-order master equation, Eq. (|13|l . and the expres- 
sion for the second-order current, Eq. (|11J) . since they are 
multiplied with exponentially small transition rates ^2L. 
As a consequence, all the needed probabilities P^\ Pj , 
are determined from Eq. (|13[) alone, which 

simplifies to 



/ V (l) V (2) 
V (l) V (2) 

4 2 T } 



\ 

Td 

(2) Jl) 
J dl ^dd 



V (2) 
y(2) 



(0) 

r 

(0) 



0. 



(15) 



plus P t (0) 



P, (0) = 1 from the normalization condition. 



T 1 ■* 1 

If one were ignorant about the described subtlety 
one might naively use the first-order master equation, 
Eq. i|12|) . with all its exponentially small (but finite) ad- 
dends to obtain a well-defined (but, in general, wrong) 
result for pP' and P) . There are situations, though, in 
which this procedure, although unjustified by construc- 
tion, leads to the correct result, namely when the total 
system is symmetric under spin reversal (nonmagnetic 
leads, p = 0), or for vanishing bias voltage, V = 0. In 
both cases, the correct result pi ^ = P^ ' = 1/2 is en- 
sured either by symmetry or as a consequence of detailed 
balance relations. It is only for broken spin symmetry 
combined with finite bias voltage V ^ that the naive 
procedure leads to wrong results. 

We remark that the current in the Coulomb-blockade 
regime far from resonance can alternatively be calculated 
without the use of the diagrammatic language. Instead 
one can employ a rate-equation approach with cotunnel- 
ing rates obtained in second-order perturbation theory 
[23l.l24l . l25j . The rate "f°',^ a for a cotunneling process, in 
which one electron leaves the dot to reservoir r' and one 
electron enters from r with the initial and final dot state 
being a and a' , respectively, is 



IT r 



^-Re J duj[l - f(uj - Mr)]/(^ - Mr') X 



pcrpcr 



pcrpcr 

f r' 



_{uj-e + i0+) 2 (u-e-U + iO+) 2 



(16) 



when the dot spin is not changed (a — a') - non-spin- fl 
cotunneling, while we get 

7rV " = ^^ Re / 'M 1 - /( W - Mr)]/(^ - Mr') X 



1 



1 



uj - e + i0+ e + U - lo + 



(17) 



for cotunneling process in which the dot spin is flipped (<r 
is the opposite spin of a) - spin-flip cotunneling. Here, 



f(uj — fi r ) is the Fermi function of reservoir r with electro- 
chemical potential fi r . The regularization +i0 + is put 
here by hand, while it naturally comes out within the 
diagrammatic formulation. There are two types of spin- 
flip cotunneling processes. Each of them involves two 
tunneling events, either through the same or through the 
two opposite tunnel barriers. Accordingly, we refer to 
them as single-barrier (r = r') and double-barrier cotun- 
neling (r r'). Double-barrier cotunneling contributes 
directly to the current, while single-barrier cotunneling 
preserves the total charge in the leads. Nevertheless, 
spin-flip single-barrier cotunneling can influence the to- 
tal current indirectly, by changing of the magnetic state 
of the dot. The probabilities P a are obtained from the 

stationary rate equation = J2rr' 
together with the normalization condition Pj + P 
The current / is, then, given by 



1. 



(7(7 



7lr 



Pa 



(18) 



This result is identical to the one obtained within the 
diagrammatic technique. Close to resonance, however, it 
is not sufficient to include the sequential and cotunneling 
processes, but also contributions associated with renor- 
malization of level position, level splitting and tunnel- 
coupling strengths become important. The diagram- 
matic language systematically takes everything into ac- 
count properly. 



C. Crossover scheme 

For both the case when sequential tunneling is allowed 
or suppressed, we have formulated a proper perturba- 
tion expansion of the current up to second order in the 
tunnel-coupling strength. When evaluating the TMR as 
a function of various parameters, such as the gate or 
transport voltage, one has to switch from one scheme 
to the other around the threshold of sequential tunnel- 
ing. At the crossover, there is no well-defined second- 
order perturbation expansion since terms of different or- 
der in r are comparable in magnitude, and their ratio 
changes continuously as a function of gate or transport 
voltage. Alternatively, we may use a crossover scheme 
that smoothly crosses over from one scheme to the other. 
This scheme consists of solving the master equation with 
first- and second-order self energies, without expanding 
the probabilities, 







'x'x 



.(2) 

J x'x 



Pv 



(19) 



and plugging this into the expression for the current, 



I =- 



2h ^ L x x 



.1(2) 

J x'x 



p. 



x ■ 



(20) 



xx' 



Up to second order in T, this result for the current is 
identical to the above-introduced accurate perturbation 
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FIG. 3: Linear conductance for nonmagnetic leads (p = 0) as 
a function of the level position. The dashed line corresponds 
to the first-order contribution G"', the dotted line represents 
the second-order conductance G*- 2 - 1 and the solid line presents 
the sum G*- 1 ' + G*- 2 - 1 . The parameters are: feT = T and 
U = 20F. The figure was generated using the scheme for the 
perturbation expansion in the presence of sequential tunnel- 
ing. 



schemes. Deviations are of third and higher order, which 
are, although unsystematic, always small for the cho- 
sen parameters, as otherwise, the perturbation expansion 
would break down anyway. 



IV. RESULTS 

A. Nonmagnetic leads 

Before presenting the results on the TMR for quan- 
tum dots attached to ferromagnetic leads, we illustrate 
the perturbation scheme introduced above for nonmag- 
netic leads. In Fig. [31 we show the linear conductance as 
a function of the level position (that can be tuned by a 
gate voltage), calculated to first (dashed line) and second 
(dotted line) order as well as the sum of both contribu- 
tions (solid line). Resonance peaks appear when either 
e or e + U crosses the Fermi energy of the leads. Away 
from resonance sequential tunneling is exponentially sup- 
pressed, and cotunneling processes dominate transport. 
But also at resonance, second-order contributions are im- 
portant, as can be seen in the figure. In particular, they 
yield a shift of the peak position and introduce an addi- 
tional broadening. 



B. Ferromagnetic leads 

We now switch to the case of ferromagnetic leads. As a 
consequence of spin-dependent densities of states in the 
leads, the dot-lead coupling strength becomes spin de- 
pendent as well. The coupling of the dot level to the 
leads acquire a factor (1 + p) or (1 — p) for coupling to 
majority or minority spins, respectively. We assume that 




FIG. 4: A sketch presenting different transport regimes. The 
respective regimes are separated by solid lines. 



spin-up (spin-down) electrons in the parallel configura- 
tion correspond to the majority (minority) electrons of 
the leads. In the antiparallel configuration, on the other 
hand, the magnetic moment of the right electrode is re- 
versed, and spin-up (spin-down) corresponds to minority 
(majority) electrons in the right lead. 

One of the main results of this paper is that the TMR 
strongly depends on the transport regime. The various 
transport regimes are sketched in Fig. 0] 

In the three diamonds around V = the number of dot 
electrons is fixed (to in regime A, 1 in regime B, and 
2 in regime A'), and sequential tunneling is suppressed. 
Sequential tunneling sets in once the bias voltage is in- 
creased above the threshold voltage, allowing for finite 
occupation of two adjacent charge states (0 and 1 for 
regime C, and 1 and 2 for regime C). In regime D all 
charge states 0,1, and 2 are possible. By performing a 
particle-hole transformation, the behavior in regime A' 
and C can be mapped to that in regime A and C, re- 
spectively. 



C. Sequential tunneling 

For reference, we list the TMR values obtained in first- 
order perturbation theory (see also Fig. |SJ. In regimes A 
(and A'), B, and D, the TMR value is 



TMR; 



A : B : D 
seq 



P 



1-p 

while for regime C (and C) it is 



2 = ^™ R ' 



lull 



TMR? eq = 



Ap 2 



P 



-TMR Ju11 



(21) 



(22) 



Within sequential tunneling the TMR through a 
quantum-dot spin valve is always smaller than Julliere's 
value for a single magnetic tunnel junction. In the lat- 
ter case, electrons are directly tunneling from one lead to 
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TMR/TMR 




eV/U 

FIG. 5: The first-order tunnel magnetoresistance as a function 
of the bias and gate voltages. The parameters are: fceT = 
1.5r, U = 40I\ and p = 0.5. 



the other. The transmission is, therefore, proportional to 
the product of the (spin-dependent) densities of states of 
both leads, i.e., proportional to (1 + p) 2 in case the spin 
of the transferred electron belongs to the majority spins 
in both leads, (1 — p) 2 in case it belongs to the minority 
spins, and p) in case it is majority spin in one 

and minority spin in the other lead. The total current for 
the parallel and antiparallel configurations is, thus, pro- 
portional to 1 + p 2 and 1 — p 2 , respectively, which yields 
Julliere's value for the TMR. 

The sequential tunneling rates in a quantum-dot spin 
valve involve the (spin-dependent) density of states of 
one lead only and are independent of the orientation 
of the other lead. To get a finite TMR, one needs to 
take into account nonequilibrium spin accumulation on 
the quantum dot, which is induced by the spin depen- 
dence of the tunneling rates. In the antiparallel configu- 
ration, the dot hosts a nonequilibrium spin accumulation 
m = (Pf — Pj.)/2 due to a different occupation of up- 
and down-spin levels in the dot, Pf =/= P[. It is, thus, 
the spin accumulation on the dot that mediates the in- 
formation about the relative magnetic orientation of the 
leads. This indirect mechanism is, however, always less 
effective than a direct coupling of the two leads, which 
is why the sequential-tunneling TMR is always smaller 
than Julliere's value. 

The result TMR = ±TMR Ju11 is characteristic of 
ferromagnet/normal-metal/ferromagnet double tunnel 
junctions without Coulomb interaction |34| . i.e., in the 
absence of any electron correlations, as well as for quan- 
tum dots with vanishing interaction C/ — > 0. For the 
regime D all three charge states play a role as for non- 
interacting case so the value of TMR also corresponds 
to this situation. The same value is reached in the 
Coulomb-blockade regimes A (A') and B, because all 



TMR/TMR' 




-1.5 -1 -0.5 0.5 1 1.5 

eV/U 

FIG. 6: The first-plus-second-order tunnel magnetoresistance 
as a function of bias and gate voltage. The parameters are 
the same as in Fig. |H| The figure was generated using the 
crossover scheme. 



transport processes in this regime are possible only due to 
hot electrons, which effectively do not feel the Coulomb 
barrier, interaction, and correlations. In regime C (C) 
Coulomb interaction is important and gives rise to the re- 
sult TMR = §TMR Ju11 . This increased TMR is related 
with the presence of a nonequilibrium spin accumulation 
and induced by it an additional charge accumulation for 
the antiparallel alignment. To illustrate this let us con- 
sider regime C for large bias voltages such that electrons 
are always entering the dot from the left and are leaving 
to the right lead. For the parallel alignment the dot oc- 
cupancy is given by Pf = Pj = Pq = h and Pd = 0, while 
the current / does not depend on the spin polarization 
p. For the antiparallel alignment, the spin-current con- 
servation condition 1° — 1^, with I" being the current 
flowing through the barrier r in the spin channel <r, yields 
(I +p)Pq = (f -p)P T and (l-p)Po = (1 +p)Pi, i.e., the 
probability Pq — (1 — p 2 ) / (3+p 2 ) to find the dot empty is 
reduced. Due to the fact that the current I ^ Pq (com- 
ing from the left lead) for both alignments, the tunnel 
magnetoresistance acquires the value |TMR Ju11 . 

As in regimes A and B sequential tunneling is exponen- 
tially suppressed, the TMR value obtained in first-order 
perturbation theory is unreliable. The TMR due to co- 
tunneling will be significantly different, as shown below. 
In regimes C and D, on the other hand, sequential tun- 
neling is present, and second-order corrections lead to 
smaller deviations only. 



D. Sequential tunneling plus cotunneling 

The TMR of first- plus second-order transport is shown 
in Fig. El where the second-order result is obtained by the 
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crossover scheme. It is clear that second-order transport 
has the strongest impact on the TMR in the Coulomb- 
blockade regime (regimes A and B). In regime B we even 
find a distinctively different behavior for the linear- and 
the nonlinear-response regimes. For regimes C and D, 
corrections due to second-order transport are smaller. 
With our theory we are able to cover all the transport 
regimes including the crossover region. In the following 
we analyze the various transport regimes in detail. 



1. Regime A 

In the Coulomb-blockade regime A the dot is empty, 
and the TMR is just due to spin-dependent non-spin-fiip 
cotunneling through the dot. There is no spin accumula- 
tion on the dot. The cotunneling rates are proportional 
to the product of the density of states of the left and right 
leads. In this regime electrons directly tunnel from one 
lead to the other similar as for a single magnetic tunnel 
junction case. Thus, the current flowing in the paral- 
lel configuration is proportional to 1 + p 2 , whereas that 
flowing in the antiparallel configuration is proportional 
to 1 — p 2 . As a consequence, the TMR is that of a single 
magnetic tunnel junction, 



2» 2 

TMR A = - „ = TMR Ju11 

1 — p z 



(23) 



i.e., twice as large as obtained within the sequential- 
tunneling approximation. 

In the regime A' the dot is occupied by two electrons 
and transport has hole- like character with only non-spin- 
flip cotunneling as for the regime A, consequently the 
tunnel magnetoresistance has the same value. 



2. Regime B 
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FIG. 7: The total linear conductance (a) in the parallel (solid 
line) and antiparallel (dashed line) configuration and the re- 
sulting tunnel magnetoresistance [solid line in (b)] as a func- 
tion of the level position. The dashed line in part (b) repre- 
sents the first-order tunnel magnetoresistance. The dotted- 
dashed curve presents the TMR calculated using the approx- 
imation Eq. Fhe parameters are UbT = 1.5r, U = 40F, 
and p = 0.5. Fhe figure was generated using the scheme for 
the perturbation expansion in the presence of sequential tun- 
neling. 



The TMR in regime B displays several nontrivial fea- 
tures. In particular, it is not constant but depends on 
both the gate and bias voltage. Furthermore, we find 
that for nonlinear response the TMR is significantly en- 
hanced as compared to linear response. In contrast, the 
TMR in the adjacent Coulomb blockade valley with even 
number of electrons, regime A, is rather trivial. This par- 
ity effect is related to the fact that the singly-occupied 
dot in regime B can be (partially) spin polarized, while 
the empty or doubly-occupied dot in regime A and A' 
respectively is nonmagnetic. 

The TMR in regime B is substantially smaller than 
that in regime A. This can be understood by the fact 
that for a singly-occupied dot both spin-flip and non- 
spin-fiip cotunneling processes are possible, in contrast 
to regime A and A' where only non-spin-flip cotunnel- 
ing occurs. There is a perfect symmetry in transmission 
magnitude between spin-flip (non-spin-fiip) processes in 
the parallel and non-spin-flip (spin-flip) in the antiparal- 
lel configuration, so in the absence of spin accumulation 



(Pf = PJ the resulting TMR would be reduced to zero. 
Only due to the presence of spin accumulation (Pj ^ P| ) 
for the antiparallel alignment transport is reduced and 
TMR > 0. Therefore, the actual value of the TMR in 
regime B depends in a sensitive way on the processes de- 
termining the spin accumulation, which is a function of 
both the gate and bias voltage. In particular, the dif- 
ferent role of spin-relaxation channels for the linear- and 
non-linear-response regime give rise to qualitatively dif- 
ferent behavior for the two cases. 

We first consider the linear-response TMR as a func- 
tion of level position (or gate voltage), as displayed in 
Fig. [7| The figure presents the linear conductance in the 
parallel and antiparallel configurations (part a) and the 
TMR (part b). We plot the first-order TMR (1) , which is 
constant and equal to half of the Julliere's value. First 
of all, one can see that the inclusion of second-order pro- 
cesses modifies the TMR substantially. The total TMR 
is well below Julliere's value as a consequence of spin-flip 
cotunneling. It is minimal in the center of the Coulomb- 
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FIG. 8: The total currents (a) in the parallel (solid line) 
and antiparallel (dashed line) magnetic configurations as a 
function of level position for eV — 20F. Part (b) shows the 
first-order contribution to the TMR (dashed line) and the to- 
tal TMR (solid line). The inset in part (b) shows the total 
TMR at lower temperature, ksT = 0.5F. The other parame- 
ters are the same as in Fig.|3 The figure was generated using 
the crossover scheme. 



blockade valley, e = —U/2, where the relative impor- 
tance of spin-flip as compared to non-spin-flip cotunnel- 
ing is strongest. To estimate the gate- voltage depen- 
dence of this relative importance we consider the ratio 
of the spin-flip over the non-spin-flip cotunncling rate, 
as given in Eqs. ((171) and iJTfij l. Since we are only inter- 
ested in the gate-voltage dependence we simply take the 
energy denominators at u> — and find that the ratio 
scales with [-l/e + l/(e + U)] 2 /[l/e 2 + l/(e + U) 2 ] = 
2/[l + (l + 2e/U) 2 }, which is maximal for e = -U/2. 
As illustrated in Fig. \7jp, the gate- voltage dependence of 
the TMR around the center is parabolic. To obtain an 
approximate analytic expression for the linear-response 
TMR, we specify our full result for the Coulomb-blockade 
regime (k-gT, T -C —e, £ + U), and take into account 
only the lowest-order corrections in the ratio x/y with 
x = \eV\, ksT, y = |e|, e + U. To describe the 
parabolic behavior, we, furthermore, expand the TMR 
up to quadratic order around e = —U/2 and obtain 



TMR 



P 



1 -p 2 



2 s 

17 



(24) 



We find that the smallest TMR value is 1/3 of that in 
regime A. As seen in Fig. [7Jd, this analytic expression 
approximates the numerical data quite well. 

We now switch to the non-linear-response regime. This 
case is illustrated in Fig. EI where the currents in the par- 
allel and antiparallel configuration as well as the resulting 
TMR are plotted as a function of the level position for 
eV — 20r. The dashed line in Fig. [HJd presents the first- 
order TMR plotted for reference. When changing the 
position of the dot level, one crosses over from regime A' 
over C to B, and then further through C to A. It can be 
seen that the behavior of TMR in regime B differs sig- 
nificantly from that in linear response, Fig. Et>. Instead 
of a minimum, we find a local maximum for e = —U/2, 
as displayed in Fig. |SJd. When lowering the tempera- 
ture, we even find a pronounced plateau of the TMR, 
with the plateau height given by Julliere's value and the 
widths determined by the region where first-order con- 
tributions are negligible. The reason for this increased 
TMR is nonequilibrium spin accumulation. The pres- 
ence of double-barrier spin-flip cotunneling, on the one 
hand, tends to decrease the TMR as discussed above. At 
the same time, on the other hand, it gives rise to spin 
accumulation that increases the TMR. As it turns out, 
the two effects compensate each other in the nonlinear- 
response regime (eV » k^,T), such that the TMR equals 
Julliere's value as if spin-flip cotunncling were absent. 
This compensation does not occur in the linear-response 
regime since in that case single-barrier spin-flip cotunnel- 
ing processes become important, which do not contribute 
to transport but reduce the spin accumulation. When 
approaching the threshold for sequential tunneling, the 
TMR drops from Julliere's value to match the first-order 
TMR' 1 '. At higher temperature, such that the plateau 
is not yet fully developed a local maximum still survives. 

The different behavior of the linear- and nonlinear- 
response regime is also nicely seen in the TMR as a 
function of transport voltage. The current for the paral- 
lel and antiparallel configuration as well as the resulting 
TMR is shown in Fig. H for e = -U/2. Unlike the first- 
order TMR' 1 ' illustrated in Fig.^D by a dashed line, the 
total TMR is a nonmonotonic function of the bias volt- 
age, which can be understood from the discussions pre- 
sented in above. For bias voltages below the threshold 
of sequential tunneling, transport is dominated by co- 
tunncling. Double-barrier spin-flip cotunneling processes 
suppress the TMR as compared to the Julliere's value. A 
finite spin accumulation, on the other hand, weakens this 
suppression and, therefore, tends to increase the TMR. 
In the linear-response regime, \eV\ <C k^T , the presence 
of single-barrier spin-flip cotunncling reduces the spin ac- 
cumulation which results in a rather low TMR. This is 
no longer the case at large bias, \eV\ 3> ksT, where only 
single-barrier spin-flip cotunncling plays no role and the 
net effect of double-barrier spin-flip cotunneling on the 
TMR is compensated. As a result we find an increase of 
the TMR in regime B with increasing bias voltage within 
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FIG. 9: The total current (a) in the parallel (solid line) 
and antiparallel (dashed line) magnetic configurations as a 
function of the bias voltage. Part (b) shows the first-order 
contribution to the TMR (dashed line) and the total TMR 
(solid line). The dotted-dashed curve presents the TMR cal- 
culated using the approximation Eq. I|26|l . The parameters 
are: k B T = 1.5T, e = -(7/2, U = 40F, and p = 0.5. The 
figure was generated using the crossover scheme. 
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FIG. 10: The differential conductance (a) for parallel and 
antiparallel configurations and the tunnel magnetoresistance 
(b) as a function of the bias voltage for different values of 
temperature. The maximum in conductance for antiparal- 
lel configuration at zero bias is clearly demonstrated. The 
other parameters are the same as in Fig. |5] Figure was gener- 
ated using the scheme for the perturbation expansion in the 
Coulomb blockade regime. 



the limits 



-TMR Ju11 < TMR B < TMR Ju11 



(25) 



The minimal value is reached at V — and e = —U/2, 
as discussed in the previous paragraph, and the maximal 
value is approached for bias voltages large as compared to 
temperature but still far away from the onset of sequen- 
tial tunneling. For an approximate analytic expression 
of the TMR around the minimum, we consider the sym- 
metric Anderson model, e — expand the TMR 
up to quadratic order in \eV\/k B T and go to the limit 
lei > k B T. The result, 



TMR B 



P 



P 



(3-pi)(eVf 



54(fc B T) 2 



(26) 



which compares well with the full numerical result, as 
can be seen in Fig. 03). When further increasing the bias 
voltage, sequential tunneling sets in. Deep in the regime 
D the TMR approaches one half of Julliere's value. As a 
consequence, the TMR has to decrease in the crossover 



regime between regimes B and D to match the correct 
asymptotic behavior, this is shown in Fig. EI 

There is one more extra feature directly at the thresh- 
old voltage for sequential tunneling. At this point, se- 
quential tunneling dominates transport but second-order 
corrections are still important. As shown in Fig. [51 this 
correction gives rise to a local minimum of the TMR as 
function of the bias voltage. To get an approximate ana- 
lytic expression for the TMR at this intersection point of 
regimes B, C and D, we assume |e| 3> k B T and expand 
the TMR up to first order in T/{k B T) to get 



TMR B|C|D = _P 

1 -p 2 

r 



i 



47rfc B T 



In 



irk B T 



, (27) 



with ^(x) being the digamma function, \&(l/2) ~ —1.96. 

The anomalous behavior of the TMR in the Coulomb- 
blockade regime is generated by the interplay of single- 
and double-barrier cotunneling for the antiparallel con- 
figuration. This is also seen in the appearance of a pro- 
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FIG. 11: The total current (a) in the parallel (solid line) 
and antiparallel (dashed line) magnetic configurations as a 
function of the bias voltage. Part (b) shows the first-order 
contribution to TMR (dashed line) and the total TMR (solid 
line). The parameters are: k B T = 1.5T, e = 20I\ U = 40I\ 
and p = 0.5. The figure was generated using the perturbation 
expansion in the presence of sequential tunneling. 



0.5 

0.4 

0.3 

S 0.2 
i— i 

0.1 
0.0 

0.65 

I 0.60 
1 

H 0.55 

^ 0.50 
H 

0.45 



0.40 




1 

■ a 


1 i 1 
parallel 




- 




* 


// 
// 


'/ \ 

antiparallel 


- 


— i — 
b 


! 1 1 1 

/ N 

' / *N. ^ 
/ / \ \ 
/ / \ * 


i — i — 
- 




' / \ * 
'/ \ 1 

\ \\ 


TMR (1) ■ 


/ — "7 
/ / 
/ / 


TMR \ 


J : 


. regime B 
■ 


regime C 
i 


regime D . 



0.5 



1.0 

eV/U 



1.5 



2.0 



FIG. 12: The total current (a) in the parallel (solid line) and 
antiparallel (dashed line) magnetic configuration as a function 
of the bias voltage. Part (b) shows the first-order contribution 
to the TMR (dashed line) and the total TMR (solid line). 
The parameters are: k B T = 1.5T, e = -lOT, U = 40r, and 
p — 0.5. The figure was generated using the perturbation 
expansion in the presence of sequential tunneling. 
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nounce zero-bias anomaly of the differential conductance 
as a function of the bias voltage in the antiparallel con- 
figuration, as we have discussed in detail in Ref. l3ll 
For completeness we repeat here some important facts 
and discuss their implications on the TMR. Deep in the 
Coulomb blockade regime such that the sequential tun- 
neling contributions can be completely ignored, we can 
use the perturbation scheme for the Coulomb blockade 
valley. In Fig. 110b we show the differential conductance 
for both the parallel and antiparallel configurations for 
different values of the temperature. For the parallel align- 
ment, the conductance shows the typical cotunneling be- 
havior, namely a smooth parabolic dependence on the 
bias voltage. This contrasts with the antiparallel con- 
figuration, for which the differential conductance has a 
pronounced zero-bias peak sitting at the bottom of a 
parabola. The width of the zero-bias peak is governed 
by temperature, indicating different spin-accumulation 
behavior for \eV\ < k B T and \eV\ > k B T. 



In Fig. ^] we show the current for the parallel and 
antiparallel configuration and the resulting TMR for the 
situation when the dot level lies above the Fermi energy 
of the leads. The first-order TMR is also shown for com- 
parison. In this case, one crosses over from regime A via 
C to D as the bias voltage is increased. At low voltage, 
regime A, current is carried by non-spin-flip cotunneling, 
with the TMR given by Julliere's value. Once the thresh- 
old to regime C is reached, sequential tunneling plays the 
dominant role. Second-order corrections to the current 
give rise to a slightly reduced TMR as compared to the 
sequential tunneling value. To find an approximate ana- 
lytic expression for this case, we consider the case of zero 
temperature, expand the TMR up to first order in T and 
assume \e\/U < 1 to get 



TMR C = 



1 -p 2 



(27 + 3 V + 3p 4 )r 



18tt(1 



p 2 )e 



(28) 



At the intersection of regimes A and C the TMR 
develops a local minimum. This is a consequence of 
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the fact that when approaching the intersection from 
regime C the sequential-tunneling-dominated TMR de- 
creases while beyond, in regime A, the TMR has to rise 
again to reach Julliere's value |35j . 

In Fig. E| we show the current as well as the first- 
order and total TMR as a function of bias voltage for 
e = — 10r. In this case, there is a crossover from regime 
B via C to D. Again, there is a local minimum of the 
TMR at the threshold to sequential tunneling due to the 
same reason as above. 



4- Regime D 

In regime D all the four dot states, i.e., x = 0, |, 4,d 
take part in transport. This situation is illustrated in Fig. 
^2f° r eV > 2(e + U). In this regime, transport is dom- 
inated by the first-order processes and the influence of 
second-order processes is negligible. Consequently, the 
value of total TMR in regime D is well described by 
Eq. J2U, as can be seen in Figs. [^b andlllb. 



E. Signature of exchange field 

It has been predicted [l9t (36| by some of us that the 
coupling of the dot levels to spin-polarized leads gives rise 
to an effective exchange field seen by the quantum dot 
electrons (an overview about the various effects of this 
exchange field is given in Ref . 13?)) . This exchange field 
is a consequence of both the Coulomb interaction on the 
dot and the spin polarization in the leads. The contribu- 
tion coming from one lead is proportional to the degree 
of spin polarization p and the tunnel- coupling strength 
r. Its direction is collinear with the leads' magnetiza- 
tion and its magnitude and even the sign is a function 
of the level position relative to the Fermi level. The to- 
tal exchange field experienced by the dot electrons is the 
(vector) sum of the two leads' contribution. This ex- 
change field gives rise to nontrivial transport behavior 
associated with a precession of the accumulated spin in 
the sequential-tunneling regime for noncollinearly mag- 
netized leads |0, 113, and leads to a split ting of the 
Kondo resonance in the strong-coupling limit |36ll38j| , as 
experimentally observed recently |14j. By applying our 
diagrammatic technique, the exchange field is automati- 
cally included. 

As we argue in the following, the exchange field will, 
under certain circumstances, also show up in the parame- 
ter regime studied in this paper, namely as an equilibrium 
spin polarization of the dot. This is distinctively differ- 
ent from the nonequilibrium spin accumulation discussed 
in the previous sections. The latter is a nonequilibrium 
effect that changes sign with bias reversal and, in partic- 
ular, vanishes for zero bias voltage. In contrast, a finite 
spin polarization at equilibrium can only occur when the 
dot level is spin split by either an external magnetic field 




FIG. 13: The occupation probabilities of the spin-up and spin- 
down dot levels as a function of the level position in the par- 
allel (a) and antiparallel (b) configuration. The zeroth-order 
occupation probabilities for the spin-up and spin-down levels 
are equal in both magnetic configurations, and are represented 
by the dotted lines. The total occupation probability of the 
spin-up (spin-down) level is presented by the solid (dashed) 
line. In the antiparallel configuration, the dashed and solid 
lines coincide. The parameters are: feT = l.ST, U = 40T, 
and p = 0.5. The figure was generated using the scheme for 
the perturbation expansion in the presence of sequential tun- 
neling. 



or by the intrinsic exchange field that we want to address 
now. 

In the antiparallel configuration, and for symmetric 
coupling to and equal spin polarization of the leads, the 
exchange-field contributions from the two leads exactly 
cancel out each other since they are of equal magnitude 
but pointing in opposite directions. This is different for 
the parallel configuration, for which the contributions 
from the two leads add up to some finite value. 

To lowest (zeroth) order in the tunnel coupling 
strengths T, the equilibrium probabilities for occupation 
with spin a =f, | are determined by the Boltzmann fac- 
tors P^ ' = P^ ' = exp(— @e)/Z, where Z denotes the 
partition function. Since the exchange field is propor- 
tional to r, it does not affect the zeroth-order occupa- 
tion probabilities, i.e., the sequential-tunneling approx- 
imation is not able to describe the exchange-field in- 
duced spin polarization. This is shown in Fig. 1131 where 
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the equilibrium probabilities calculated to zeroth- and 
zcroth- plus first-order in the dot-lead coupling are pre- 
sented. A finite spin polarization for the parallel config- 
uration is only generated by the first-order corrections 
7^ Pi ■> that we obtain by solving the master equa- 
tion given by Eq. I|13|) . The e-dependence of the spin 
polarization seen in Fig. 1131 reflects the e-dependence of 
the exchange field. The exchange field for a particle-hole 
symmetric band vanishes in the middle of the Coulomb 
blockade valley, e — —U/2, and has different sign on ei- 
ther side. As a consequence the dot polarization changes 
sign as well. 

Since in regime B S^^ are exponentially suppressed, 

the exchange splitting and probabilities P^ do not affect 
the second-order transport. These probabilities affect 
only higher-order transport contributions, which at low 
temperature T <T K lead to the Kondo effect HHMUl. 



V. SUMMARY 



We have discussed electronic transport through quan- 
tum dots coupled to ferromagnetic leads. Based on a 
formalism that allows for a systematic perturbation ex- 
pansion in the tunnel coupling strength, we analyzed the 
TMR through a single-level quantum dot for the linear- 
and nonlinear-response regime, at or off resonance, with 
an even or odd dot electron number. We found different 
TMR values for different transport regimes. In addition 
to the full numerical results we provided approximate 
analytic expressions for various limiting cases. The most 
important findings are: 

(i) Except for the Coulomb-blockade valley with an even 
dot-electron number and the nonlinear-response regime 
of the Coulomb-blockade valley with an odd dot-electron 
number, the TMR is below that of a single magnetic tun- 
nel junction. 

(ii) There is an even-odd asymmetry between the 
Coulomb-blockade valleys with an even or odd number 
of electrons, that is related to the absence or presence of 
spin-flip cotunncling, respectively. 

(iii) In the Coulomb-blockade valley with an odd number 
of electrons, the TMR values for the linear and nonlinear 
response regimes differ strongly from each other, associ- 
ated with different spin-relaxation processes that affect 
the spin accumulation. 

(iv) The linear-response TMR in the Coulomb-blockade 
valley with an odd number of electrons is a function of 
gate voltage, which reflects the relative importance of 
spin-flip and non-spin-flip cotunncling. 

(v) The TMR at the onset of sequential tunneling dis- 
plays a local minimum, which is a consequence of inter- 
polating the TMR behavior away from resonance. 
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APPENDIX A: DIAGRAMMATIC TECHNIQUE 

In this Appendix we present general rules in energy 
space for calculating contributions of various diagrams. 
We also present an exemplary calculation of one of the 
second-order self-energies. Afterwards, we show how to 
determine self-energies contributing to electric current. 

1. Rules in energy space 

Contribution of a particular diagram to the self-energy 
S x / X can be found following the general rules in the en- 
ergy space: 

1. Draw all topologically different diagrams with fixed 
time ordering and position of vertices. Connect the 
vertices by tunneling lines. Assign the energies of 
respective quantum dot states to the forward and 
backward propagators. To each tunneling line as- 
sign a frequency w, the spin of tunneling electron 
and label of the junction. 

2. Tunneling lines acquire arrows indicating whether 
an electron leaves or enters the dot. For tunnel- 
ing lines going forward with respect to the Keldysh 
contour assign a factor 7,7 <t (cj), whereas for tunnel- 
ing lines going backward assign 7 J t <T (w). 

3. For each time interval on the real axis limited by 
two adjacent vertices draw a vertical line inside the 
interval and assign a resolvent 1/(AE + i0 + ), with 
AE being the difference of all energies crossing the 
vertical line from right minus all energies crossing 
the vertical line from left. 

4. Each diagram gets a prefactor (— l) b+c , with b be- 
ing the number of vertices lying on the backward 
propagator and c denoting the number of crossings 
of the tunneling lines. 

5. Each internal vertex represents a matrix element 
(xl-^lx')i w ith A being a dot operator, A = d^,d a . 
Consequently, a minus sign may appear due to 
these matrix elements. This is because |d) = 
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d£ = — dt | cr) (depending on the definition of 
state |d)), where tr =f or a = L To account for this 
factor, multiply each diagram by (— l) m , where m is 
the number of vertices connecting the spin-cr state 
with doubly occupied state. 

6. Integrate over all frequencies and sum up over the 
reservoirs. 

The parameters 7^ t<T (w) are defined as 

= S-^-^ (Ai) 

7 r _ » = ?[l-/(^-Mr)], (A2) 



with f(x) being the Fermi-Dirac distribution function, 
f(x) — 1/ [ex.p(x / ksT) + 1], and [i r representing the 
electrochemical potential of lead r. 



2. Calculation of 

In order to find the zeroth-order and first-order prob- 
abilities, one needs to determine all the self-energies of 
first and second order in V. Below, we present an exem- 
plary calculation of one of the second-order self-energies, 

(2) (2) 

. The equation for Yr s „ can be graphically presented 

as 




To calculate the self-energy, it is necessary to evaluate ing the general rules described above, the corresponding 
each contributing diagram. As an example, we present contribution, A3, is given by 
calculation of the third diagram of Eq. I|A3(1 . Follow- 



A, 



E 



dui du 2 %° {lo\ )7+ <t (u>2 ) 



1 



U + i0+ L02 



iO- 



<A4) 



The first (second) factor on the right-hand side follows 
from the rule 4 (5). There are also three resolvents ac- 
cording to the rule (3). Among the various diagrams 

(2) 

contributing to S^, there is a diagram (eleventh in 
Eq. IA3(1 whose contribution is equal to minus complex 
conjugate of the contribution due to the third diagram, 
An = — Rc(A.3) + Tm(A3). This can be shown by in- 
terchanging the backward and forward propagators and 
changing the direction of the tunneling lines. As a conse- 
quence, the real parts of these diagrams cancel, whereas 
the imaginary parts add to each other. Thus, it is nec- 



essary to determine only the imaginary part of one of 
those two diagrams, A3 + An = 2Tm(A3). After contour 
integration, the imaginary part of A3 is given by 



Im(A 3 ) = £ £ 



per 

--p^h^Ti + Hr 2 - £a ~ £g - U)X^ r a 2 (2fl r2 -£g-U) 



r CT 



(A5) 
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with /b (x) being the Bose-Einstein distribution function 
Ib(x) = 1/ [exp(x/k-QT) — 1]. The corresponding coef- 
ficients A^{e a ) are defined as, A^(e a ) = X^(e a ) — 
X±?{e a + U), with X±°{e a ) = ±T° /{2it)B a {e a - M r) 
and B a {x) given by 



B a+1 (x) 



Re 



27rfc B T 



-In 



W 



27rfc B T 



where ^(z) is the digamma function, and we have used 
the Lorentzian cutoff function of the form p v {^>) = 



W 2 /[(w — /i„) 2 + W 2 ], with W being the cutoff parame- 
ter. As contribution from a single diagram may depend 
on W, the final result does not. In the calculations the 
cutoff parameter was taken to be equal to 100r. 



In a similar way, one can calculate contributions of all 
diagrams, which give 



.(2) 

~* (7 <J 



-2m \l^{e a )Xir^) 

ri,r 2 I 



)+ 7 +^)X 2 --( £ff ) 



+7,7 (e„ + U)X+l(e s + U)+ j+f(e s + U)X^{e a + U) 



Z7T 



X£ (e a ) + X+l {e a +U) + -A+°(e 



r CT 

' 2tt 



(A6) 



3. Diagrams contributing to the current 

To find current flowing through the system, one has to 
determine the self-energies S 1 , see Eq. <|10|) or (|llf) . This 
can be done by realizing that each term of the expansion 
of the current operator / is equal to the corresponding 
expansion term of the reduced density matrix multiplied 
by a factor of e/H. The only difference is that now for 
each external vertex lying on the upper (lower) branch 
of the Keldysh contour, corresponding to tunneling of 
an electron into the left (right) or out of the right (left) 
lead, we have a multiplicative factor +1/2, whereas for 
each external vertex on the upper (lower) branch of the 



contour, describing tunneling of an electron into the right 
(left) or out of the left (right) lead, there is a factor of 
-1/2. 

We have determined all the first-order and second- 
order self-energies contributing to electrical current, S 1 ' 1 ^ 
and S 1 ^, and found that from the first-order self- 

i vJW ^(ri y^W • 

energies only E y, £ CT \/, £ CT y, £ d y give nonzero con- 
tributions. In the case of the second-order self-energies 

1(2) 

we found S xx = 0, with x — 0, T>l>d. This is how- 
ever only the case for the current operator defined as 
I = (/r — Il)/2, where I r is the current operator for 
electrons tunneling to the lead r. 
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